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ABSTRACT 

With galaxy groups constructed from the Sloan Digital Sky Survey (SDSS), we analyze 
the expected galaxy-galaxy lensing signals around satellite galaxies residing in different 
host haloes and located at different halo-centric distances. We use Markov Chain 
Monte Carlo (MCMC) method to explore the potential constraints on the mass and 
density profile of subhaloes associated with satellite galaxies from SDSS-like surveys 
and surveys similar to the Large Synoptic Survey Telescope (LSST) . Our results show 
that for SDSS-like surveys, we can only set a loose constraint on the mean mass of 
subhaloes. With LSST-like surveys, however, both the mean mass and the density 
profile of subhaloes can be well constrained. 

Key words: cosmology: dark matter - galaxies: haloes - methods: statistical - galax- 
ies: subhalo - gravitational lensing 



1 INTRODUCTION 

In the cold dark matter (CDM) scenario, large-scale struc- 
tures in the universe grow hierarchically through gravita- 
tional instabilities. Galaxies are assumed to form in dark 
matt er potential wells through gas cooling and star forma- 
tion l|White fc Reedl 19781 : IWhite fc Frenklll99ll ). During the 
hierarchical formation process, when small haloes merge into 
larger systems, they become subhaloes. High resolution sim- 
ulations show that while some of them are disrupted dure 
to processes such as tidal stripping and impulsive heating, a 
large fraction of the subhaloes survive. Hence, probing the 
masses and density profiles of the population of subhaloes is 
a key test for the CDM structure formation model. 

The mass function, spatial distribution and density pro- 
file of subhaloes has been extensively studied with semi- 
analytical models a s well as A^-body simulations (e.g., 
Havashi et al.l l2003l: [Cao et all 1 2004': 'van don Bos ch et all 



On the other hand, it is very challenging to probe 
dark matter subhaloes observationally because of their dark- 
ness and the relatively weak gravitational potential com- 
pared to that of their host haloes. Arguably the best (and 
most direct) probe of dark matter substructure is gravita- 
tional lensing. The existence of substructure in a smooth 
dark matter halo induces flux-ra tio anomalies for mul- 
tiple images of a lensed sys tem (Ma o fc Schneider 19981: 
Metcalf fc Madaull200ll:lMao et al...2004l:lKochanek fc Daial 



2004 iMaccio fc Mirandal l2006l : IXu et all |2009| ). and also 



20051: iGiocoH et all l200i 



2010; Springcl ct al 



2008 



Zentner fc Bullock 20031: iTavlor fc Babui2004l:IOguri fc Led 



2004 iGill et al. ,2004') State-of-the-art, high-resolut ion sim- 
ulations ( Springel et al.l [20081 : iDiemand et al.ll2007l ) can re- 
solve subhaloes down to a mass of ~ lO'^h-^Me, and thus 
provide detailed predictions for both their mass function and 
density profiles. 
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perturbs the surface brightness of extended Einstein rings 
and arcs (Ko opmans 2005; Ve getti fc K oopmans 2009a|b|; 
IVegetti et abi |201G| . |2012| ). So far, about 200 galaxy- 
sized strong lensi ng systems have been discovered (e.g., 
iBohon et ;il]l2008t ). The constraints on the mass fraction of 
subhaloes in galaxies have been investigated (e.g.. IXu et al.l 
[iooab However, the number of strong lensing systems with 
high quality imaging observations is still limited. Further- 
more, strong lensing effects can only probe the very central 
region of galaxies. Therefore it is not easy to obtain a general 
understanding about subhaloes from strong lensing effects 
alone. 

Because most of the satellite galaxies are thought to 
reside in subhaloes, galaxy-galaxy lensing can be an effec- 
tive way to probe subhaloes statistically. While it was first 
developed to estimate the dark matter distribution of mas- 
sive systems, the recent advance of wide and deep surveys. 
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such as the Sloan Digital Sky Survey (SDSS) and the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) 
0, has allowed applications of galaxy-galaxy lensing analy- 
ses to study the mass distributions around lens galaxies of 
different luminosi ties, stellar masses, c olors, and morpho- 



logical types (e . g..lBrainerd et al.ll 19961 : iHudson et al] 
Hoekstra et all l2003l: iHoekstral |2004| : iMcKav et al " 



Mandelbaum et al.'"'2005', 2006", 200^; ISheldon et ~ 



1998 



2001 



200£ 



Johnston et al.. ,2007 '). Several studies have investigated 



the potential of using galaxy-galaxy lensing to probe the 
masses and dens i ty profiles of dark matter subhaloes(e. 
lYang et al.l l2006l : iLi et all l2009l : iPastor Mira et all [201 
Current observations can only set partial constraints on 
subhalo properties for individual massive clusters of galax - 
ies (e.g.. (Limousin et al.|[2007l : iNataraian et al.|[2007l . l2009l ). 
However, with the next generation of large surveys, such 
as the Large Synoptic Survey Telescope (LSST) 0, the sur- 
face number density of source galaxies that can be used for 
galaxy-galaxy lensing analyses can reach Ug ~ 50 arcmin"^, 
in comparison with ~ 1 arcmin"^ for SDSS. This will sig- 
nificantly increase the signal-to-noise of the lensing signal of 
dark matter subhaloes, thus enabling direct measurements 
of their masses and density profile. 

The goal of this paper is to examine the potential of us- 
ing galaxy-galaxy lensing to constrain the properties of dark 
matter subhaloes, such as their mass and density profile. 
The subhalo properties of satellite galaxies with some fixed 
property (i.e., stellar mass) are likely to depend on both host 
halo mass and the location of the satellite galaxy within that 
host halo. In order to be able to probe these dependencies, 
we need to distinguish satellite galaxies located in different 
haloes and at different radii from the centers of their host 
haloes. One way to do t his is to select lens galaxies using 
a galaxy group ca talog (' Yang et al. I l2006l : IJohnston et al.l 



I2OO7I : ISheldon et a l. 2009; Li et al.l bOoi T This allows one 



to select as lenses a subset of satellite galaxies that reside 
in haloes (groups) of similar mass, and that are at simi- 
lar (projected) distances from their halo (group) center. In 
iLi et al.l (|2009l l , we applied such a method to predict galaxy- 
galaxy lensing effects for lens galaxies of different luminosi- 
ties and differe nt morphological ty pes using a group catalog 
constructed by lYang et al.l l|2007l . hereafter Y07) from the 
SDSS. The predictions are fou nd to agree well wi t h lens - 
ing observations of SDSS from [Mandelbaum et al.l l|2006h . 
demonstrating the validity of the method. In this paper, 
we use the same methodology to predict the galaxy-galaxy 
lensing effects for satellite galaxies selected from the SDSS- 
DR7 group catalog. We investigate the corresponding sig- 
nal detectability with current and next generation surveys. 
Employing a Markov Chain Monte Carlo (MCMC) method, 
we further explore the possibility of constraining both the 
subhalo and host halo properties in lensing observations at 
different noise levels expected from different surveys. 

This paper is organized as follows. We provide a brief 
description of the galaxy-galaxy lensing basics in Section [21 
and discuss the modeling method in Section [S] In Section 
13.11 we introduce the group catalog, SDSSGC, from which 
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lens galaxies are selected. In Sections l3.2l andl 3.3l we describe 
our models for the dark matter distribution around galaxies. 
In Sections |3] and [SI we show the results and examine the 
detectability of the predicted lensing signals in SDSS-like 
and LSST-like surveys. We discuss some systematic bias in 
our method, and how to correct for it, in Section [S] Section 
[7] contains a summary. 

Throughout the paper, we adopt a ACDM cosmol- 
ogy with parameter s given by the WMAP-7-year data 
([Komatsu et al.ll2010l '). 



2 GALAXY-GALAXY LENSING 

Galaxy-galaxy lensing measures the tangential shear, 7t(i?), 
azimuthally averaged over a thin annulus at the projected 
radius R around the lens galaxies. In the weak lensing 
regime, this quantity is related to the excess surface den- 
sity, AE (hereafter ESD) through the relation 



AS(7?) = 7t(i?)Ecrit = E(< R) - E{R) , 



(1) 



where E(< R) is the average surface mass density within 
R, and E(J?) is the azimuthally averaged surface density at 
R. It is noted that there is no mass-sheet degeneracy here, 
and AE(i?) is independent of a uniform background. In the 
above equation. 



Ecrit — 



47rG DiDisil + zi) 



(2) 



is the critical surface density in comoving units, with Ds 
and Di the angular diameter distances to the lens and to 
the source, Dis the angular diameter distance between the 
lens and the source, and zi the redshift of the lens. 

The lensing signal around a galaxy is determined by the 
projected density profile around it. On average, the surface 
mass density, E(_R), is related to the line-of-sight projection 
of the galaxy-matter cross-correlation function, 5g,m('"). Un- 
der the approximation that lenses are at distances much 
larger than R, we can write 



and 



E(< R) = ^ J T.{u)udu, 



dx ; 



(3) 



(4) 



where p is the mean density of the universe and x is the 
comoving radial distance along the line of sight. 

The lensing signal around a satellite galaxy depend 
sensitively on its location in t he host dark matter halo 
(jYang et al.ll2006l : iLi et al.1 120091 '). The E(i?) around a cen- 
tral galaxy, which mostly resides at the center of the host 
dark matter halo, is dominated by the density profile of its 
host halo. On the other hand, the lensing signal of a satel- 
lite galaxy, which orbits around in the host halo, consist of 
two parts. On small scales, the signal is dominated by the 
subhalo associated with the satellite itself. On larger scales, 
however, the lensing signal is mainly due to the host halo. 
We therefore need to model the density profiles of both host 
haloes and subhaloes. 

When calculating the surface mass density around 



3 



a satellite, we neglect the contributions from other sub- 
haloes. This approximation should not lead to large er- 
rors because the fraction of mass contained in sub- 
haloes is only about 10% of the to t al mass of th e host 



halo (e.g. van den Bosch et all l2005l : ISpringel et al] l2008l : 



iGiocoli et al.ll201(]| ). For a single halo, these substructures 
produce small fluctuations on the host halo profile. However, 
in galaxy-galaxy lensing analysis, one stacks the lensing sig- 
nal around many lens galaxies, and the net contribution from 
subhaloes other than the ones associated with the satellite 
galaxies themselves is averaged out and included in the host 
halo profile. In the model calculation, we also neglect the 
two-halo term, specifically, the contribution to the lensing 
signal from other halo es in the foreground and background . 
Our previous studies l|Li et al.l l2009l : ICacciato et"al] 120091 ') 
calculate this contribution with different methods, and show 
that the two-halo term is completely negligible on the scales 
of individual haloes we are concerned with here. 



3 MODELING THE STRUCTURE OF DARK 
MATTER HALOES 

In this paper, we adopt the same methodology as that used 
m iLi et al.l (|2009l ) to model the galaxy-galaxy lensing signal 
around a sample of satellite galaxies. In the following sub- 
sections, we describe briefly the galaxy and group catalogs, 
and our models for the mass distributions around satellite 
galaxies. 



a halo mass to a group is that there is a one-to-one relation 
between Mateiiar and the halo mass. Once a theoretical dark 
matter halo mass function is adopted, one can establish a 
relation between halo mass, M, and Mstoiiar so that the num- 
ber of haloes with masses above M is equal to the number 
of groups with characteristic stellar mass above Mstoiiar- A 
group with a given Msteiiar is then assigned the correspond- 
ing halo mass M. Clearly, this one-to-one mapping requires 
the group sample to be complete. Therefore, we only use 
complete samples of groups in the SDSSGC in our ranking. 
The masses of other groups are estimated using linear in- 
terpolation based on the Mstoiiar -M rela tion obtained from 
the c omplete sample. We refer readers to I Yang et al.l (|2007l . 
l2008h for details about the group catalog construction and 
the halo mass assignment. According to Y07, the uncertainty 
in the mass assignment is about 0.2-0.3 dex for groups con- 
sidered in this paper. This uncertainty will not change our 
results significantly because the considered lensing signals 
are the average signals over a statistical sample of galaxies. 
We have tested the effect by performing calculations with 
group mass to which an artificial 0.3 dex log-normal error 
is added. This uncertainty brings negligible change in our 
results. Note that the halo mass assigned to a group in the 
SDSSGC is M200, which is the mass enclosed in the radius, 
r2oo, defined such that M200 ~ 47rr2oo(200p)/3. For consis- 
tency, we convert M200 to our definition of halo mass (see 
Eq. [7] belo w) using the conversion method described in the 
appendix of IHu fc Kravtsovl l|2003l '). 



3.1 Galaxy groups 

In iLi et all ||2009|'). w e used the SDSS DR4 group cata- 
log IjYang et al. I HootI '). Here we use an updated version of 
this catalog B (hereafte r SDSSGC) based on the SDSS DR7 
l|Abazaiian et al.l I2OO9I ). The group catalog is constructed 
with the a daptive halo-based group finder developed by 
lYang et al.l (^005, 2007) using galaxies with spectroscopic 
redshifts in the range of 0.02 ^ 2: ^ 0.2. The redshift com- 
pleteness is C > 0.7. Three group samples with different 
sources of galaxy redshifts have been constructed. Our anal- 
ysis is based on Sample II, which consists of 599301 galaxies 
with redshift from the SDSS and 3269 galaxies with red- 
shift from other sources. There are in total 472113 groups, 
including those with only one member galaxy. 

A key aspect of this group finder is to estimate the 
halo mass, M, for each group with a ranking method. In 
SDSSGC, two estimators for halo mass are provided. One 
is based on the characteristic luminosity of a group, de- 
fined to be the total luminosity of all member galaxies with 
Mr — 5 log ft < 19.5. The other is based on the characteris- 
tic stellar mass, Mstoiiar, defined to be the total stellar mass 
of members galaxies with Mr — 5 log ft < 19.5. The stellar 
mass of an individual galaxy is calculated from it s luminos- 
ity an d colors using the fitting formula given et al.l 
l|2003h . Y07 showed that the characteristic stellar mass is a 
better indicator of the halo mass, and thus we adopt this 
mass estimator throughout the paper. 

The basic assumption of the ranking method to assign 



http: / /gax. shao.ac.cn/data/Group. html 



3.2 Host halo density profile 

We assume that the host dark matter halo of each group is 
centered on the most massive group member, to which we 
refer as the central galaxy. T he dark matter host haloes are 
assumed to follow the NFW (|Navarro et al.lll997l ) profile. 







a X ^ c 

else 



where x = r/vs, with being the characteristic scale of the 
halo and related to the halo virial radius r-vir through the 
concentration parameter c = r-vir/'"s, and pcrit is the critical 
density of the universe. The characteristic over-density So is 
related to the average over-density of a virialized halo, Avir, 
by 



So 



Avir 



3 ln(l+c)-c/(l-hc) 



(5) 



For the ACDM model considered here, we adopt the para- 
metric form of Avir given by Bryan & Norman (1998) based 
on the spherical collapse model. 



Avir = IStt^ + 82[Q{z) - 1] - 39[Q{z) - 1]^ 



(6) 



where Q{z) is the cosmological density parameter at redshift 
z. The viral mass of a halo can then be written as: 



47r 3 

M = — rvirAvirPcrit ■ 



(7) 



It is seen that given the mass of a halo, its density profile 
depends only on the concentration parameter c. Numerical 
simulations show that at a g iven redshift, c decr e ases grad- 
ually with halo mass (e.g., iBullock et all I2OOII : lEke et all 
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However the exact mass dependence of the concen- 
tration parameter has not yet been well constrained by 
observations. Various fitting formulae have been proposed 
on the bases of numerical simulations (e.g., Bullock et alj 



200ll:IZhao et al ll2003l : iDolag et al.l2004l : iMaccio et al.ll2007l : 
Zhao et al.ll2009l ). Since the difference between these differ- 
ent fitting functions does not af fect our results quali tatively, 
we adopt the fitting formula of iBulIock et"ai] (|200ll ). 



M 



1 + z V lOi^h-iMg 



(8) 



where c, ~ 8 for the ACDM cosmology considered here. 
Although simulations indicate a s catter of ^ O.ldex in this 
concentration mass relatio n (e.g. Jind 2000l: Bullock et al.l 
I2OOII : IWechsler et al. I I2OO2I : iMaccio et al.ll2007l ). we ignore 
this scatter in our analysis. We have verified that adding 
scatter has no significant impact on any of our results. 

To obtain galaxy-galaxy lensing signals around a halo 
center, we need to projec t the 3-D mass distribution. Ac- 
cording to iHamana et al.l l|2004l ). the ESD can be written 
as 



AE(y) = 2AE./(j/),j/ = 



R 



(9) 



where AEs = paVs- For the NFW profile, the dimensionless 
function f(y) can be written as (| Wright fc Brainerd|[2OO0l : 
iBartelmannfliggel ) 



fiy) = { 



1 - 



iiy<l 

if y = 1, 
if y > 1 



(10) 



Note that the above equation assumes that the NFW profile 
extends infinitely without a truncation at the virial radius. 
We have tested that the difference between truncated and 
non-truncated profiles has negligible impact on our results. 

3.3 Subhalo density profile and mass distribution 

In order to assign a subhalo mass to a s atellite galaxy we use 
the ev olved subhalo mass function of Ivan den Bosch et al.l 
which describes the abundance of subhaloes as func- 
tion of their evolved, present-day mass. Specifically, for each 
satellite galaxy, we assume that the mass of its subhalo at 
accretion (i.e., its un-evolved mass) is related monotonically 
to its stellar mass. The retained mass fraction of the sub- 
halo after evol ution in the host halo can be described by a 
parameter fm ■ iGao et all (|2004 ) analyzed the radial depen- 
dence of fm from a large sample of subhaloes in cosmological 
simulations, and showed that 



/m = 0.65(rdis/rvir)^''^ : 



(11) 



where rdis is the distance of the subhalo from the center 
of the host halo, and rvir is the virial radius of the host 
halo. Observationally, the line-of-sight distance cannot be es- 
timated accurately, and only the projected halo-centric dis- 
tance, rp, can be used. We then adopt the following approach 
to estimate the three-dimensional distance of the satellite. 
For a satellite with a given rp, we randomly sample a 3-D 
halo-centric distance assuming that the spatial distribution 



of the satellites follows the NFW profile, and use this dis- 
tance in equation to estimate fm- With fm, we can 
define a ranking parameter Q for every member satellites in 
a group with host halo mass Al as 



Q = fmM, 



(12) 



where A/* is the stellar mass of the satellite galaxy. We then 
generate a set of subhalo masses for a given host halo mass 
M using t he fitting formula for the ev olved subhalo mass 
function of Ivan den Bosch et al.l (|2005l ). Finally, by ranked 
mapping, a subhalo mass is assigned to a satellite according 
to its ranking parameter Q. 

FigH] shows the results of our subhalo mass assignment. 
The upper panels are the stellar mass distributions for the 
satellites with their host halo mass in the range indicated in 
the plots. We split satellites into sub-samples according to 
their projected halo-centric distance rp shown by different 
lines. One can find that the peak of the stellar mass distribu- 
tion does not depend strongly on their projected distance. 
The mass distributions of the correspondng subhaloes are 
shown in the lower panels. The average mass to light ratio 
is about 10. 

We model the density profile of subhaloes with a trun- 
cated NFW profile 



Psub{r) 



ftPi,sub{r) if r ^ rt, 
if r > Vt 



(13) 



where Pi,sub{r) is the NFW profile corresponding to the mass 
of the subhalo at the time of its accretion into its host. 
The parameter ft is a dimensionless factor describing the 
reduction in the central density, and rt is a cut-off radius 
imposed by the tidal force of the host halo. In Pi.sub{r), 
the characteristic scale and density are denoted by rs,sub 
and So^sub, respectively. Note that the parameters of ft and 
So, sub can be combined to a single parameter denoted by 
/!>o,sub. For ft = 1 and rt » rs,sub, psubir) approaches to the 
standard NFW profile pi,sub{r)- For the cut-off radius rt , we 
use the analytical tidal ra dius formula (|Binnev fc Tremaind 
ll987l : lTormen et aLliTooi ). 



n 



Me, 



(2 - d In M/d In r)M(< rdis 



1/3 



rdis , 



(14) 



where M{< rjis) is the host halo mas s with in a sphere of 
radius rdis. As shown bv lSpringel et af] (|2008l ). this analyti- 
cal prediction agrees well with the trunctation radii of dark 
matter subhaloes A''-body simulations. The density profile is 
normalized to the mass assigned to the subhalo by choosing 
a proper ft (or equivalently po,sub) . Therefore in our model 
the mass profile assigned to a subhalo is specified by three 
quantities: (i) the stellar mass of the satellite galaxy; (ii) the 
host halo mass; and (iii) the distance between the satellite 
and the center of the host halo. 

It should be pointed out that there are still substan- 
tial uncertainties in modeling the mass distribution around 
individual satellite galaxies. In particular, many of the re- 
sults about subhaloes are obtained from pure A^-body sim- 
ulations. It is unclear how significant the effect of including 
baryonic matter is. However, the aim of this work is to in- 
vestigate to what extent current and future lensing data can 
constrain the density distribution of dark matter subhaloes. 
Our relatively simple model for the density distribution of 
subhaloes should be sufficient for this purpose. 
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Figure 1. The figure shows the stellar mass and subhalo mass distribution of satellite galaxies with host halo mass in certain ranges. 
The host haloes mass range are [10^'^, 5 X lO^'']h~^M0 for left panels and [10^**, 5 X 10^^]h~^MQ for right panels. In the upper panels 
of the figure, we plot the stellar mass distribution for the satellites. We split the satellites into sub-samples according to their projected 
halo-centric distance rp and plot the stellar mass distribution of these sub-sample with different color. The respective rp range for each 
color of the line is marked in bottom left panel in unit of r^i^. In the lower panels, we plot the mass distribution for the subhalocs for 
corresponding satellites sample. 



4 MODELLED LENSING SIGNAL 
4.1 Lensing signal of individual satellite 

We first calculate the behavior of galaxy-galaxy lensing sig- 
nal as a function of the projected radius arou nd a single 
satelli te. Si milar res u lts ca n be found in, e.g. lYang et al.l 
l|2006t) and iLi et al.l (|2009l ). In FiglH we plot AE(_R) for 
satellites of different mass and at different position in a 
host halo with mass of M = lO^'^h'^M©. In the left panel, 
the subhalo mass is set to be zero to show the lensing sig- 
nal contributed by the host halo alone. Different lines show 
the predictions for satellites located at different projected 



halo-centric distances, Vp. It is clear that the contribution 
from the host halo depends strongly on the position of 
the satellite. For Vp — 0, i.e., the central galaxy, E(7?) is 
just the projection of NEW density profile, and AS(7?) de- 
creases monotonously. For a satellite galaxy with rp 7^ 0, 
however, AE(_R) from the host halo is nearly on small 
scales around the satellite. This is because the host halo 
density varies smoothly on small scale around the satellite. 
As R grows, AE(7?) decreases to negative values, reaching 
a minimum at J? = rp where the outer annulus in Eq.((T)) 
reaches the centre of the host halo. It then goes up rapidly, 
eventually approaching the AE(_R) profile for the central 
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galaxy. In the right panel, we show AE(_R) for satellites at 
Tp = 0.5h~^Mpc with different subhalo mass. Clearly sub- 
halo dominates the inner part of the ESD profile. The value 
of AE(i? = 0.01h~^Mpc) increases by a factor of 2.5 when 
the subhalo mass increases from 

4.2 Stacking 

Observationally, weak lensing signals are derived by accu- 
rately measuring the shape of the light distribution of source 
galaxies. The intrinsic shape of source galaxies then con- 
tributes significant noise to the lensing signal. Specifically, 
the measured tangential ellipticity e+ of a source galaxy is 
related to the lensing tangential shear 74 acting upon it by 

e+ = 27t7^ + eV"^ (15) 

where e™* is the intrinsic tangential ellipticity of the source 
galaxy and TZ is the "responsivity" , refiecting how the 
shape of a galaxy resp onds to the shear applied to it 
(IBernstein &: Jarvi 3 I2OO2I '). This quantity can be determined 
from observational data, and we set TZ = 0.87 following 
iMandelbaum et al.l (|2006l ). In galaxy-galaxy lensing, one 
suppresses the noise arising from the intrinsic shape of 
source galaxies by combining shape information from as 
many source images as possible. Unless the density of back- 
ground source is extremely high, this typically requires also 
stacking the signal from many lens galaxies. Assuming no in- 
trinsic alignment for source galaxies, the averaged intrinsic 
ellipticity over many galaxies is expected to approach zero. 
Thus the average of e+ gives rise to an unbiased estimate of 
2'ytlZ. The corresponding uncertainty of the tangential shear 
measurement can be written as 

2na^ = ^/^l~{~7i/y^N~, (16) 

where Npair is the number of lens-source pairs, asN ~ 0.3 
is the source galaxy intrinsic shape dispersion for one com- 
ponent of the ellipticity, and ae is the measurement noise 
for one component of the ellipticity. The measurement noise 
can originate from photon noise and inaccurate PSF cor- 
rections. For the SDSS, falls in the range from 0.05 
to 0.4, depending on the lum inosity of the source galax- 
ies (|Mandelbaum et al.ll2005l ). Throughout this paper, we 
adopt (Te — 0.2. It should be pointed out that the LSST 
will improve significantly in survey depth and angular res- 
olution compared to SDSS, and thus is expected to be 
much smaller. The measurement noise adopted here is there- 
fore a very conservative estimate for future surveys. We fix 
the lens redshift to be zi = 0.15, which is the mean redshift 
of galaxies in the SDSSGC. For source galaxies, we con- 
sider two models. The first model (hereafter LEVI) is for 
SDSS-like surveys, which, to a certain extent, represents the 
current state-of-the-art for large surveys. The second model 
(hereafter LEV2) is for future LSST-like surveys. For LEVI 
and LEV2 the source galaxy redshift is taken to be Zs = 0.3 
and Zs — 1, respectively. For simplicity, we do not consider 
detailed redshift distributions for source and lens galaxies. 
This may lead to inaccurate predictions for lensing signals, 
especially for LEVI with relatively low Zs- However, for the 
purpose of comparing the detectibillity of LEVI and LEV2 
surveys, our simplification should be adequate. On the other 
hand, for future studies requiring high precision, the redshift 



distributions of source and lens galaxies has to be properly 
accounted for. 

With fixed zi and Zs, we have 

crAE(-R) = a^,{R) X T,crit{zi,Zs) . (17) 

Thus the lensing measurement noise only depends on A'^pair, 
which in turn is determined by the number of lens galax- 
ies and the number density of source galaxies. We use the 
SDSSGC to estimate the number of lens galaxies. FigO 
shows the number of lens galaxies in SDSSGC as a func- 
tion of the projected group-centric radius for different halo 
mass range. Typically the SDSSGC provides between 2000 
and 6000 satellite galaxies {— lenses) per bin in host halo 
mass and group-centric radius, for the binning adopted here. 
When halo mass decreases, the number of groups keeps in- 
creasing, but the number of satellites per host goes down. At 
halo mass range of 10^^ — lO^*h~^M0 , one obtains largest 
number of lens galaxies. 

Future surveys such as LSST will include spectroscopic 
data, making it difficult to construct a reliable group cat- 
alog from the survey data itself. This will have to await 
future deep and wide spectroscopic surveys which will allow 
the construction of very large group catalogs, out to high 
redshifts. This would allow the galaxy-galaxy lensing based 
subhalo studies proposed here to be extended to higher red- 
shifts, using tomography. An alternative is to consider us- 
ing photometric redshifts for the construction of group cat- 
alogs. We leave these avenues for future investigation. In 
what follows we simply adopt the lens statistics from the 
SDSSGC for both LEVI and LEV2; i.e., even for LEV2, we 
still only consider groups constructed from relatively shal- 
low SDSS-like spectroscopic surveys. Finally, for the number 
density of source galaxies we adopt 1.6 arcmin"'^ for L EVI 
(as appropriate for SDSS, see IMandelbaum et al.|[2005l ) and 
60 arcmin"^ for LEV2. 

In Fig|4l we compare the lensing measurement noise for 
LEVI and LEV2, respectively. We calculate average lensing 
signals around satellite galaxies with host halo mass in the 
range of [10",5 x 10^*]h-^Mo and [10^^ 5 x 10^=^]h-^Mo 
and with halo-centric distance in the range of [0.1,0.2] and 
[0.5, 0.6] h~^Mpc. We plot the expected noise from our two 
noise models, LEVI and LEV2, with red and blue error 
bars, respectively. It is seen that both the SDSS-like survey 
and the LSST-like survey can detect the lensing signal from 
the host halo well. For the inner parts where the subhaloes 
dominate, the observational noise from SDSS-like survey is 
much larger than the signal. On the other hand, a LSST-like 
survey can detect the signal with high S/N. The S/N does 
not drop for smaller groups, because the number of smaller 
groups is much larger than that of massive ones and the 
mean subhalo mass does not drop significantly in smaller 
groups (see Fig[l]). 



5 MODEL INFERENCE WITH THE MARKOV 
CHAIN MONTE CARLO METHOD 

5.1 Markov Chain Monte Carlo fitting 

For a given set of observational data 9, the posterior prob- 
ability of the model parameters vr can be derived from the 
likelihood function, C{6\tt), and the prior probability, -P(vr), 
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Figure 2. In left panel, we plot the lensing signal around satellites at different positions in the halo of lO^^h^^MQ. The satellite mass 
is set to zero. Different lines represent different projected halo-centric distance. In the right panel, different lines represent lensing signal 
around satellites of different mass within a host halo of 10-'^^h~^MQ. 
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Figure 3. The figure shows the number of galaxies in SDSSGC as function of projected group-centric radius. Different line styles 
represent different host halo mass. The respective ranges of log(M/h~^M0) are marked in the figure. 
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Figure 4. In this figure, we plot the galaxy-galaxy lensing signal around satellites in certain host halo mass bins and projected distance 
bins and compare it with the noise estimation. The red errorbars show the noise estimation of LEV1(SDSS like survey), while the blue 
rectangles show the LEV2 noise(LSST like survey). 



of these parameters. According to Bayes' rule, we can write 



Pie) 



(18) 



where the normahzation, P(0), is called the evidence. The 
prior probability P(7r) describes our known knowledge about 
the model parameters. In our fiducial computation, we sim- 
ply adopt a fiat prior over a range in the parameter space. 
As a test of the robustness of our inferences, we also use 
another set of priors. We write the likelihood function C as 



ln£ : 



E 



\ CTAE 



(19) 



where AE(_Ri) is the observed excess surface density at ra- 
dius bin Ri and AE(_Ri|7r) the theoretical prediction with 
model parameters tt. In this paper, we use MCMC to ex- 
plore the posterior distribution P{tt\6). The key component 
of the MCMC method is a guided random walk in the pa- 
rameter space. The likelihood function at a certain volume of 
the parameter space is then proportional to the number den- 
sity of points in that volume. The MCMC sa mpler used here 
is that provided in the CosmoMC package (|Lewis &: Bridle! 



|2002|). which adopts the Metropolis-Hastings algorithm 
(iMetropolis fc UlamI 1 19491 ; [Me tropolis ot al. 1953; Hastingi 
Il97d ) b y default. A detailed revie w of this method can be 
found in IChib fc Greenberd l| 19951 ). For each fitting process, 
we generate three MCMC chains starting from different po- 
sitions in the parameter space. We use the runtime conver- 
gence criteria in CosmoMC, which computes the standard 
Gelman and Rubin i?-statistic diagnostics to monitor the 
convergence. We declare convergence when R < 1.1. Only 
the second half of the chain values are used to sample the 
posterior probability. 

5.2 Model inference 

Here we investigate to what extent the observations of satel- 
lite galaxy-galaxy lensing from LEVI and LEV2 surveys can 
constrain the satellite and host halo properties described by 
a set of parameters, including (M, c, Vp) for the host halo and 
(Afsub, Po,sub, T's.sub) for the subhalo. We construct the 'ob- 
served' galaxy-galaxy lensing signals for satellites selected 
from SDSSGC following the descriptions in §4. We consider 
separately two subsamples of satellites according to their 
host halo mass, one in the range of [10^'', 5 x lO^'']h~'^M0 
(SAMPLEl), and the other for [10^^ 5 x 10^^]h-^MQ (SAM- 
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Table 1. The mean values of parameters for the input galaxies, 
and the boundary of the parameter space adopted in the MCMC. 
M and Mgut, are in units of h~^M0; rp and sub are in units of 
h-^Mpc; Po.sub is in units of IO^^/i^MqMpc-^. 





logM 


c 




log A4ub 


PO.sub 


^s,sub 


mean input 


14.247 


6.89 


0.55 


11.67 


0.99 


0.0155 


high bound 


14.5 


10 


0.6 


12 


10 


0.2 


low bound 


13.5 


3 


0.5 


9 


0.1 


0.001 



PLE2). The projected halo-centric distances of satellites are 
chosen to be in the ranee of [0.5, 0.6]h~^Mpc for both sub- 
samples. We stack the galaxy-galaxy lensing signals for the 
satellite galaxies in each subsample thus constructing a set 
of 'observed' signals. We then employ MCMC fitting to these 
mock data in order to derive constraints on the mean values 
for the six parameters, three for host haloes (M, c, r^) and 
three for subhaloes {Ms^b, po.sub, ''s.sub)- The selection of the 
bin size in the projected halo-centric distance of satellites is 
a compromise between statistical errors and systematic bias. 
If we used the bin size of 0.2h~^Mpc instead of 0.1h~^Mpc, 
for example, the source galaxy number would increase by 
about a factor of two and thus the statistical errors are de- 
creased. However, that stacking satellites over a larger range 
in rp leads to larger biases in the derived mean host halo 
mass M and concentration parameter c from MCMC fitting 
to the stacked satellite galaxy-galaxy lensing signals. 

To derive the parameters we are interested in, we fit the 
mock lensing data with the model. 



T.{R) = Ehost(-R|M,c,rp) + Ssub(iilAfsub,Po,sub,rs,sub) 



(20) 



where 



Ehost(i?|M, c,rp) 



1 

2^ 



E(^r2 + i?2 + 2TpRcose)de, 
(21) 

where E(J?) is the projected density profile of the NFW 
host halo with mass M and concentration parameter c. Note 
that R is the distance to the satellite galaxies around which 
we detect galaxy-galaxy lensing signals. The second term 
Ssub(-R) is the projected density of the subhalo with profile: 



Psub 



if r < rt 



(l + r/r,_,ub)^('-/'-=,Bub) " ■ " ■ ' ^ (22) 

if r > rt 



where the truncation radius, rt, is computed using Eq. (|14p . 

We first analyse SAMPLEl with relatively massive host 
haloes. Table [1] shows the mean values of the six parameters 
of the input sample, and the boundary of the fiat priors we 
adopt in our MCMC fitting. 

We separately analyse the constraints expected from 
LEVI surveys and LEV2 surveys. FigOshows the marginal- 
ized posterior probability distribution for the model parame- 
ters for LEVI. The contours indicate the 68% and 95% confi- 
dence levels. The last panel in each row shows the marginal- 
ized probability distribution for the corresponding parame- 
ter, in which the probability distribution (blue histogram) 
and the average value (vertical line) from the input sample 
are also shown for comparison. As is evident, even for LEVI 
the host halo mass M and concentration parameter c can 



already be constrained reasonably well. However, the con- 
straints on the subhalo mass are extremely weak, with the 
68% confidence range covering two orders of magnitude. For 
the subhalo density profile, no meaningful constraints can 
be obtained from LEVl-type surveys. 

The results obtained using the LEV2 noise model are 
shown in Fig[6] Clearly, the factor 50 increase in the number 
density of source images causes a dramatic improvement in 
the constraints on the model parameters compared to those 
in FigO The host halo mass and the concentration param- 
eter can be constrained with high precision. The subhalo 
mass can also be tightly constrained with a la confidence 
range of about 0.2 dex. For the amplitude and scale radius 
of the subhalo density profile, reasonable constraints can be 
achieved. Note, though, that there is a strong degeneracy be- 
tween po,sub and rs,sub, as seen from their joint constraints 
in the fifth panel of the bottom row. 

To see more clearly the improvements from LEVI to 
LEV2, in Fig[7]we directly compare the marginalized prob- 
ability distributions for M, c, Msub and rs,sub obtained us- 
ing the LEV2 (solid) and LEVI (dashed) noise models. This 
demonstrates the remarkable potential of the next genera- 
tion of LSST-like surveys for subhalo studies using satellite 
galaxy-galaxy lensing. 

We also perform the MCMC fitting for SAMPLE2 with 
group sized host haloes with mass in the range of [10^'' , 5 x 
10^^] h-^M0. The results are shown in Fig[S] It is seen that 
for group sized host haloes and the satellites therein, we can 
still get good constraints with LEV2 surveys. 



5.3 The impact of prior choice 

The results above are obtained with fiat priors for all pa- 
rameters. Here we test the impact of prior choice on our 
inferences of model parameters. Since the groups used here 
are selected from SDSSGC group catalog, each group has 
already been assigned an estimated mass. Thus we do have 
some idea about the mass distribution of the selected groups 
that may be used as priors in the MCMC fitting. As an 
approximation, we model the mass distribution with a log- 
normal function, 



P(log M) 



2-Ka 



, (logM- log M)^ \ 



where M is the mean mass of the selected groups, and 
o" ~ (o"o + o\)l\fN , with (70 = 0.3 the mass assignment 
uncertainty in the SDSSGC, cri the standard deviation of 
the group mass in logarithmic space, and A'^ the number of 
selected groups. In FigEl we show the constraints using this 
distribution as the prior for the host halo mass. The left pan- 
els show the marginalized posterior distribution of M, and 
the right panels are for the subhalo mass Msub- The upper 
and lower panels are for LEV2 and LEVI, respectively. The 
black and blue histograms are the results with flat priors for 
all the parameters, and the lognormal prior for M and flat 
priors for the other parameters, respectively. The red lines 
in the left panels illustrate the lognormal prior distribution 
for M . For LEV2, because the constraints are already tight, 
adding the lognormal prior on M does not change the con- 
straints signiflcantly, although it does decrease the bias in 
M by a small (barely significant) amount. In the case of 
the LEVI noise model, the posterior distribution for M is 
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Figure 5. The marginalized posterior probability distribution for the model parameters in LEVI case. The contours show 68% and 95% 
confidence levels. The last panel of each row shows the 1-d marginalized probability distribution, together with the average value of the 
input satellites (vertical solid lines). The blue histograms show the corresponding distributions from the input SAMPLEl. 



essentially identical to its prior distribution, indicating that 
the LEVI lensing data does not improve the constraints on 
host halo mass beyond our prior knowledge. For the subhalo 
mass, the posteriors based on both LEVI and LEV2 are not 
affected by the prior on M. This is expected from Figs[5]and 
[6l which show that the degeneracies between host halo mass 
and subhalo parameters are very weak. 



6 BIAS AND CONTAMINATIONS 
6.1 Stacking bias 

From Figs[5]and[6l it can be seen that the peak of the pos- 
terior distribution for some parameters deviates from the 
input mean value. These biases are from fitting the model 
to galaxy-galaxy lensing data obtained by stacking a large 
sample of satellite galaxies.. The properties of the satellite 
galaxies and their host haloes in the sample are not identi- 
cal but spread over certain ranges. The parameters obtained 
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Figure 6. The same as FigObut for results using LEV2 noise model. 



from the fitting to tfie stacked signals then correspond to the 
results of certain averages over such a sample. Depending on 
the quantities, the averages obtained may be biased relative 
to the means of the input values. In our analyses, although 
the properties of lens galaxies, such as their host halo mass 
and the halo-centric locations, are selected to be similar, 
they still cover finite ranges with some distributions. In par- 
ticular the subhalo mass of the lens galaxies covers a very 
broad range. Bias arises simply because the lensing signal 
depends on these parameters in a non-linear fashion. For 
subhalo mass, the difference between the peak of the poste- 
rior distribution and the input mean value is about 0.1 dex. 



slightly larger than the 2a width of the posterior distribu- 
tion. 

One way to suppress the bias is to use narrower mass 
bins. In our model, the subhalo mass of a satellite depends 
mainly on its stellar mass although the host halo mass and 
halo-centric radius also influence somewhat through tidal in- 
teractions. Too see the effects of different binning, we split 
satellites in our fiducial satellite sample into five finer stel- 
lar mass bins. The corresponding lensing signal of each of 
the mass bins is shown in Fig llOl respectively. We then can 
analyze the host halo and subhalo properties with MCMC 
by fitting to the lensing signals from the five subsamples 
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Figure 7. A close comparison between tlie posterior distributions shown in Figs[5] and |6] The solid histograms are the marginalized 
distribution of M, c, Mgub and gut, for LEV2, while the dashed histograms for LEVI case. 



jointly. In principle, we should consider five sets of parame- 
ters (Mgub, rs_sub, Po.sub) each for a single bin. We then need 
to deal with 15 parameters for subhaloes, plus the ones for 
host haloes. The task would be difficult with the statistics of 
the expected data. On the the hand, one naturally expects 
certain relations between different quantities, which can be 
parameterized with a much smaller number of free parame- 
ters. By MCMC fitting to the 'observed' data, we can extract 
the constraints on these parameters. This approach allows 
us to control the number of free parameters and at the same 
time to model the lensing signals better than that from a 



single broad bin of subhalo mass. As a test, we assume that 
Msuh/M, depends on the stellar mass M, according to a 
power-law, 

Me,b f M, y 
-Mr^°°l lO"h-^Me J ■ (24) 

Furthermore, we assume that Vs^suh and rt depend on sub- 
halo mass through the following parameterizations 
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Figure 8. The same as Fig[7l except that here results are shown for sateUites in host haloes with masses in the range [10^^, 5x 10^^]h ^Mq 
and with projected halo-centric distance in range [0.5, 0.6]h~^Mpc. 



and 



n = Co + ci n^suh , 



(26) 



where ao, ai, ^'o, 6i, cq and ci are all free parameters. 
We then fit the mock lensing data from all 5 stellar mass 
bins simultaneously to derive constraints on these param- 
eters. Fig llll shows the 68% and 95% confidence range of 
the marginalized posterior probability distributions for these 
subhalo parameters. For comparison, the 'true input' values, 
obtained by directly fitting relations ([24]), (|25|) . and (|26p to 
the input subhaloes, are marked with the plus symbol in each 



panel. We can see that the 68% posterior contours enclose 
the input values, indicating that the bias due to binning can 
be effectively reduced by dividing lens galaxies into fine bins. 



6.2 Contamination from fake group members 

So far we have not considered possible contaminations in the 
group catalog and assumed that all galaxies assigned to a 
group are true members. In reality, however, some galaxies 
that are identified as satellites may actually be central galax- 
ies of other (low-mass) haloes along the line of sight. In what 
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Figure 9. The figure shows the effeet of priors of host halo mass on parameter constraining. The black(blue) histogram represents the 
marginalized posterior distribution of M(left column) and subhalo mass(right column) with flat prior (lognormal prior). The upper panels 
show the results for LEV2 case and the lower panels show the results for LEVI case. In the left panels, the red solid lines show the 
prior distribution of halo mass. We scale the amplitude of prior distribution, so that it can be shown in the same panel with posterior 
distribution. 



follows, we refer to such galaxies as interlopers. For galax- 
ies of the same luminosity, central galaxies produce much 
stronger lensing signals than that of satellites. Hence, even 
an interloper fraction of 10% can introduce significant errors 
in the inferred subhalo parameters. 

To estimate the impact of interlopers, w e make use 
of the SDSS mock group catalog provided by lYang et ahl 
This mock catalog is constructed by running the 
halo-based group finder of lYang et al. I (|2005l ) on a mock 
galaxy redshift survey (MGRS) built by populating dark 
matter haloes with galax ies according to the cond itional 
luminosity function(CLF; Ivan den Bosch et al.1 120071 1. The 
CLF, which describes the luminosity distribution for galax- 



ies in haloes of a given mass, is constrained using the clus- 
tering and abundances of galaxies in the SDSS. As the re- 
sult, the luminosity function and clustering properties of the 
MGRS accurately matches those of the SDSS. . The MGRS 
also takes into account of real observational conditions by 
mimicking the s ky coverage and co mpleteness trend of the 
SDSS survey ("see lYang et al.l |20o3) for the detail). For such 
a mock group catalog, we not only know the group to which 
a galaxy is assigned, but also the dark matter halo to which 
the galaxy truly belongs. Thus, it is particular suitable to 
examine the potential bias that arises due to interlopers. 

From the SDSS mock group catalog, we select satellites 
in groups with assigned host mass from the group finder in 
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Figure 10. The lensing signal around satellite galaxies in different stellar mass range. All satellite galaxies are selected from groups 
of [10^^,5 X lO^*]h~^M0, and with projected halo centric distance of [0.5, 0.6]h~^Mpc. Each subplot shows lensing signal of satellites 
in certain stellar mass bin. We mark the log(Af * /h~^MQ) range for each subsample in subplots. The red errorbars show the noise 
estimation of LEVI (SDSS like survey), while the blue rectangles show the LEV2 noise (LSST like survey). 



the range of [10^*, 2 x lO"]h"iM0 and with the projected 
halo-centric distance of [0.5, 0.6]h~^Mpc . We use the model 
described in Section [3] to generate mock galaxy-galaxy lens- 
ing data. To isolate the errors due to interlopers, we fix the 
subhalo mass to be lO^^'^h~^M0 in calculating the fiducial 
data, which is similar to the mean subhalo mass of SDSSGC 
satellites used in previous sections. We find that about 10% 
of the selected satellite galaxies are actually centrals of other 
host haloes (i.e., are interlopers). We then assign each of 
them a NFW profile according to their host halo mass. The 
resulting lensing signals are shown in Fig |12l where the solid 
line is the excess surface density profile of the true satellite 
galaxies, and the dashed line shows the signal due to the 
interlopers. As can be seen, the interlopers contribute about 
15% of the total signal in the inner part, which in turn can 
lead to large bias in the model fitting. Fig |13l shows the bias 
in subhalo mass introduced by these interlopers. We fit two 
sets of lensing signals separately. The first set contains con- 
tributions from both true satellites and interlopers. This is 
referred to as the 'mixed' signal. The second set contains 
only true satellite galaxies, and is referred to as the 'true' 
signal. The solid line is the constraint from the 'true' signals 
and the dashed line is from the 'mixed' signals. The vertical 
line indicates the input subhalo mass. We can see that the 
subhalo mass inferred from the 'mixed' signal is biased high 
by ~ 50% compared to that inferred from the 'true' signal. 
The latter result is very consistent with the input value. 

The above analyses shows that it is important to care- 
fully consider the impact of interlopers in the group cata- 



log. Unfortunately, it is virtually impossible to completely 
eliminate interlop ers. In fact, the halo based group finder of 
lYang et al.l (|2005l l has been optimized to minimize the im- 
pact of interlopers, among some other constraints. Rather 
than trying to reduce (or eliminate) interlopers, one may 
also try to account for them in the modeling. Using empiri- 
cal relations between the luminosity /stell ar mass of centra l 
galaxies and their host halo mass (see e.g. lYang et al.ll2007l ) 
it is fairly straightforward to fit for subhalo mass and in- 
terloper fraction simultaneously. As a test, we assume that 
each galaxy in the selected sample has the same possibility 
Peon to be a central galaxy. We use the relation between the 
central galaxy luminosity and its host halo mass from our 
MGRS to assign a halo mass to a central galaxy. The total 
signal is then modeled as. 



AE(ii) = AEsat(-R)(l - Peon) + ASccn (ii)P« 



(27) 



where AEsat is the lensing signal calculated by assuming 
no contamination, and ASccn is the lensing signal calcu- 
lated by assuming that all satellites selected are actually 
central galaxies of other haloes. The fitting result is shown 
in Fig |14l Although the posterior distribution shows degen- 
eracy between Peon and Mgub, the bias due to contamination 
is now suppressed. If the MGRS is sufficiently realistic, we 
should in principle be able to obtain some estimates for the 
interloper fractions, which can then be used as a prior in 
the MCMC fitting. This should allow for tight and unbiased 
constraints on the subhalo mass. 
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Figure 11. Fitting results for lensing signal in Fig llOl Contours show 68% and 95% confidence range of tlie marginalized posterior 
probability distribution. The plus marks show the true inputs value. 



7 SUMMARY 

In this paper, we have investigated the feasibihty of con- 
straining the masses and density profiles of dark matter 
subhaloes associated with satellite galaxies using galaxy- 
galaxy lensing. With the use of a group catalog constructed 
from a large redshift survey, such as the SDSS, we can ef- 
fectively distinguish central and satellite galaxies. Therefore 
we can select satellite galaxies according to their host halo 
mass, halo-centric distance and stellar mass. By stacking the 
galaxy-galaxy lensing signal of satellite galaxies with simi- 



lar properties, we can then study both the host halo and 
subhalo density profiles. In this paper, we have used the 
SDSS group catalog constructed by Yang et al. (2007) to 
predict the galaxy-galaxy lensing signal for satellite galax- 
ies residing in different host haloes and located at different 
halo-centric distances. To examine to what extent such data 
can be used to infer the properties of host and sub-halo, we 
have considered two different noise levels, LEVI and LEV2, 
that correspond to the levels of measurement noise expected 
from a current generation galaxy survey, such as SDSS, and 
from a next generation galaxy survey, such as LSST, respec- 
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Figure 12. The figure shows the lensing signal of galaxies 
in mock SDSS group catalog. The galaxies are selected from 
haloes of [10i'*,2 x lO"]h-iM0 with halo-centric radius of 
[0.5, 0.6]h^^Mpc. The signal from true satellites is represented 
by solid line and that from fake members by dashed line. 
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Figure 13. The 1-d constraint on the subhalo mass. The dashed 

and solid lines are for the results derived from fitting to the 
'mixed' and 'ture' signals of satellites, respectively. The vertical 
line indicates the input value of the subhalo mass. 



Figure 14. Fitting to the mixed lensing signal of mock catalog. 
The contour shows 2d posterior distribution of Mguj, and Peon- 
The plus mark represents input value. 



tively. Using the MCMC method, we investigated how well 

LEVI and LEV2-type data can constrain the mass and den- 
sity profile for subhaloes and host haloes simultaneously. For 
satellite galaxies in massive groups, with host halo masses 
in the range of [10^*, 5 X 10^'*]h-iMo, the density profile of 
the host halo can be well constrained for both LEVI and 
LEV2 noise levels. However, the data quality that is achiev- 
able with current LEVl-type surveys is insufficient to put 
any meaningful constraints on the subhalo properties. In the 
case of a LEV2-type survey, on the other hand, the galaxy- 
galaxy lensing data is predicted to be of sufficient quality 
that one can put tight constraints on the average subhalo 
mass, with a la confidence of about 0.2 dex. Even the ampli- 
tude and characteristic scale of the subhalo density profiles 
can be constrained, albeit with still relatively large uncer- 
tainties. We also demonstrate that, with LEV2-type surveys, 
it is even possible to probe subhaloes in group-sized host 
haloes with masses as low as lO"h-^M0. 

We also discussed some potential systematics that re- 
sult in biased estimates. One of these arises from the fact 
that one stacks the signal from satellite galaxies that span 
a significant range in properties of host halo and subhalo. 
Since the lensing signal does not scale linearly with model 
parameters, the best fit of the mean value of the parameters 
can be biased relative to the underlying values of the stacked 
sample. We have shown that such bias can be reduced by 
dividing the satellite sample into finer stellar mass bins and 
using parameterized forms for the scaling relations between 
satellite and subhalo properties. Another bias arises from 
the presence of interlopers in the group catalog (i.e., from 
galaxies that have erroneously been assigned to a group due 
to projection effects). This implies that some of the galaxies 
identified as satellites in the group catalog are actually cen- 
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trals of other (typically low-mass) haloes. Our test using a 
mock SDSS group catalog shows that about 10% of the satel- 
lites are such interlopers. Such a contamination can bias the 
inferred subhalo mass high by ~ 50%. To overcome the bias 
effect, wc propose to include the contamination fraction as a 
free parameter in the model fitting. Our test shows that the 
bias in the subhalo mass due to the contamination can be 
effectively removed at the expense of enlarged uncertainties. 
This uncertainty, in turn, can be reduced by using priors on 
the interloper fractions which can be obtained from realistic 
mock galaxy redshift surveys. We conclude that measure- 
ment of galaxy-galaxy lensing around satellite galaxies with 
future surveys such as LSST holds great promise for con- 
straining the properties of dark matter substructure. 
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